Spontaneous Transition of Turbulent Flames to Detonations in Unconfined Media 
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Deflagration-to-detonation transition (DDT) can occur in environments ranging from experimen- 
tal and industrial systems to astrophysical thermonuclear (type la) supernovae explosions. Substan- 
tial progress has been made in explaining the nature of DDT in confined systems with walls, internal 
obstacles, or pre-existing shocks. It remains unclear, however, whether DDT can occur in unconflned 
media. Here we use direct numerical simulations (DNS) to show that for high enough turbulent 
intensities unconfined, subsonic, premixed, turbulent flames are inherently unstable to DDT. The 
associated mechanism, based on the nonsteady evolution of flames faster than the Chapman- Jouguet 
deflagrations, is qualitatively different from the traditionally suggested spontaneous reaction wave 
model, and thus does not require the formation of distributed flames. Critical turbulent flame 
speeds, predicted by this mechanism for the onset of DDT, are in agreement with DNS results. 

PACS numbers: 47.70.Pq, 47.40.Rs, 97.60.Bw 



Since the discovery of detonations, the question of the 
physical mechanisms that create these self-supporting, 
supersonic, shock-driven reaction waves has been a fore- 
front topic in combustion theory. Uncontrolled develop- 
ment of detonations poses significant threats to chemical 
storage and processing facilities, mining operations, etc. 
[H, while controlled detonation initiation in propulsion 
systems could revolutionize transportation [2|. On as- 
trophysical scales, detonation formation is presently the 
most important, yet least understood, aspect of the ex- 
plosion j^, 01 powering type la supernovae, which, as 
standard cosmological distance indicators, led to the dis- 
covery of the accelerating expansion of the Universe @] . 

Early studies 0] showed that a detonation can arise 
from a slow, highly subsonic defiagration ignited in an 
initially unpressurized system. Significant progress has 
since been made experimentally Q and numerically 0- 
in elucidating the physics of the deflagration-to- 
detonation transition (DDT) in confined systems, and 
particularly in closed channels. These studies showed 
that the confining effect of channel walls on the hot, ex- 
panding products of burning and the interaction of the 
resulting fiow with walls and obstacles are important in 
accelerating the fiame and causing the pressure increase, 
thus creating conditions necessary for the detonation ig- 
nition. This raises the question: Is DDT possible in un- 
confined media without assistance of walls or obstacles, 
e.g., in unconfined clouds of fuel vapor or in the interior 
of a white dwarf star during a supernova explosion? 

Zel'dovich et al. [isj originally suggested that a deto- 
nation can form in a region ( "hot spot" ) with a suitable 
gradient of reactivity. The resulting spontaneous reac- 
tion wave propagating through that gradient creates a 
pressure wave that can eventually develop into a shock 
and then a detonation [l3.[l5l|. In confined systems, mul- 
tidimensional direct numerical simulations (DNS) have 
shown that hot spots can form through repeated shock- 



fiame interactions and fuel compression by shocks [lo| . 

It remains unclear, however, if and how hot spots could 
form in unconfined, unpressurized media. The most 
likely mechanism involves fiame interactions with intense 
turbulence. It was sugg ested [Hill] that the flame struc- 
ture could be disrupted by turbulence, producing a dis- 
tributed fiame with reactivity gradients capable of initi- 
ating a detonation. There are, however, no realistic ah 
initio experimental or numerical demonstrations of this 
process. Here we show that high-speed turbulence-flame 
interactions can indeed lead to DDT, but through a dif- 
ferent process in which pressure build-up in the system 
does not rely on the propagation of global spontaneous 
reaction waves and, thus, does not require the formation 
of large-scale gradients of reactivity. 

Model and method. — The DNS presented here solve 
the compressible reactive-flow equations including ther- 
mal conduction, molecular species diffusion, and energy 
release IT], |l8| . They use an ideal-gas equation of state 
and a single-step, first-order Arrhenius kinetics to de- 
scribe chemical reactions converting fuel into product. 
Simplified reaction-diffusion models represent stoichio- 
metric Il2-air and CH4-air mixtures with unity Lewis 
number and reproduce both experimental laminar fiame 
and detonation properties 
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Simulations were per- 
formed with the code Athena- RFX, which uses a fully un- 
split corner-transport upwind scheme with PPM spatial 
reconstruction and the HLLC Riemann solv er [17,-, 2ftl. 
Turbulence is driven using a spectral method |17|, |2l| . 

Numerical simulations. — Fig ure 1 shows a traditional 
combustion regime diagram [22| with a summary of the 
cases studied. Regions of the diagram representing dif- 
ferent burning regimes are bounded by lines of constant 
nondimensional Damkohler, Da, Karlovitz, Ka, and 



Reynolds, Re, numbers [22] ■ Cases 6, 7, and 10 represent 
several simulations testing numerical issues such as the 
solution convergence and the absence of unphysical ef- 
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FIG. 1. (Color online) Combustion regime diagram [24] show- 
ing the simulations discussed here. Symbol color and shape 
indicate the reactive mixture and the mode of burning. The 
full flame width /f,o ~ 2(5l,o [13. 



fects due to the boundary conditions. All simulations are 
well-resolved with the resolution at least Ax — (5l.o/16, 
where 5l,o cm is the thermal width of the laminar flame 
in cold fuel [l7, 18|. Convergence was confirmed for Cases 
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18[ and 7 using Ax — (5l,o/8 — (5l_o/32, and con- 
vergence during the DDT process was confirmed in Case 
10 for Ax = (5l,o/8 - (5l,o/16. 

This paper focuses on Case 10 and later compares it to 
other simulations shown in Fig. 1. Case 10 is a DNS of 
a premixed H2-air flame interacting with the high-speed, 
steadily driven turbulence. Its setup is similar to the pre- 

, which analyzed a 
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vious detailed study of Case 6 
steady turbulent flame evolution in a smaller system with 
lower intensity turbulence. The computational domain is 
a uniform 256 x 256 x 4096 Cartesian mesh with width 
L — 0.518 cm, giving the resolution Ax — (5l^o/16 with 
^Lfl ~ 0.032. Kinetic energy is injected at the scale L 
to produce homogeneous, isotropic turbulence with the 
characteristic velocity U — 1.^ x lO'' cm/s ~ 63S'l^o at 
the scale L, where Sl^o — 3.02 x 10^ cm/s is the laminar 
flame speed in cold fuel. The large-scale eddy turnover 
time is = L/U = 27.3 /iS, the integral velocity is 
Ui = 1.2 X 10** cm/s fvAOSLfi, and the integral scale is 
Z = 0.12 cm. Resulting turbulence away from the flame 
has an equilibrium Kolmogorov energy spectrum ocfc~^/^ 
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in the inertial range extending to scales < Sl,o 

Initially, fuel has temperature Tq = 293 K and pres- 
sure Pq — 1.01 X 10^ erg/cm^. Steady-state turbulence 
is allowed to develop for 2Ted- At this point {t — 0) 
a planar flame is initialized normal to the x-axis. The 
boundary conditions are zero-order extrapolations at the 
x-boundaries and periodic conditions at the y- and z- 
boundaries. After « 2Ted, the turbulent flame is fully 
developed and reaches a quasi-steady state (QSS) that 
lasts until t w 6.5Ted- Figure [5] shows the turbulent flame 
speed, St, based on the fuel-consumption rate [l7| . Tur- 



FIG. 2. (Color online) The turbulent flame speed, St, nor- 
malized by the sound speed in cold fuel, Ca,o- The legend gives 
the Case numbers (Fig.[l]). The two shaded gray regions show 
the range of critical values of St (see eq. U])) based on the 
sound speed in fuel, c^j, and product, Cs.p, for fuel temper- 
atures in the range 360 — 430 K. Blue dots on the curve for 
Case 10 indicate times of individual profiles in Fig. [S] Time 
is normalized by the corresponding value of T^d in each case. 



bulent flame properties during this period are consistent 
with the earlier analysis of such QSS in Case 6 [l7i [is}. 
In particular, the flame folded inside the flame brush re- 
mains in the thin reaction-zone regime with its reaction 
zone structure virtually unaffected by turbulence. St is 
primarily controlled by the increase of the flame surface 
area with an additional periodic increase < 30 — 40% due 
to flame collisions and the formation of cusps. 

In contrast to Case 6, the QSS in Case 10 is relatively 
brief (Fig. After t w 6.5ret;, St begins to increase 
rapidly, becoming supersonic by 7.18Ted and exceeding 
the Chapman- Jouguet (CJ) detonation velocity, Dcj, at 
7.5Ted- DDT occurs at 7.53Ted, and St reaches its maxi- 
mum at 7.58Ted- At 7.63rerf, a fully developed overdriven 
detonation emerges and quickly relaxes to Dcj- 

The system evolution during this process is shown in 
Fig. El At 6.39Ted, a slight overpressure arises inside the 
flame brush, but the energy-generation rate per unit vol- 
ume, E, is still close to its value in the planar laminar 
flame. As the pressure grows and the turbulent flame ac- 
celerates, fuel inside the flame brush is compressed and 
heated. This increases the local flame speed, Sl, causing 
E to rise. At later times, E exceeds the laminar value 
by 2 orders of magnitude. Such accelerated burning 
leads to further fuel compression and larger Sl- The re- 
sulting feedback loop drives a catastrophic runaway pro- 
cess that produces a large pressure build-up and creates 
strong shocks inside the flame brush. These, in turn, cre- 
ate conditions in which a detonation can arise. (Details 
of this last stage will be presented in a separate paper.) 

Up until the moment of DDT, the average fuel tem- 
perature, Tf, inside the flame brush remains < 700 K 



3 




FIG. 3. (Color online) The j/-z-averaged profiles of (a) pres- 
sure, P, (b) fuel mass fraction, Y, (c) x-velocity, Ux, (d) 
energy-generation rate per unit volume, E, and (e) tempera- 
ture of pure fuel {Y > 0.95) in Case 10. The time from igni- 
tion for each profile is shown in panel (b) and indicated with 
blue dots in Fig. [2] i? is normalized by its value in a planar 
laminar flame propagating in cold fuel, El,o = iSl,oPo/5l,o, 
where q is the chemical energy release and po is the density 
of cold fuel. T[ is shown inside the flame brush up to the 
moment of detonation formation (t < 7.58red). 



(Fig. [3^), and the corresponding induction times are 
much larger than all dynamical timescales. At all times, 
the average internal flame structure (reconstructed using 
method described in 17[) is close to that of a laminar 
flame in fuel with the corresponding Tf and pressure. 
Thus, during the runaway, burning is controlled by flame 
propagation and not by autoignition, which precludes the 
formation of global spontaneous reaction waves. 

Mechanism of the spontaneous runaway. — Consider 



an unconfined fluid volume V with the total internal en- 
ergy e. To increase the pressure inside V (as in Fig. [3]), 
an energetic process must generate energy comparable 
to e on the characteristic sound-crossing time of this vol- 
ume, i.e., i ~ s/ts- If this volume represents a flame with 
width 5 and cross-sectional area L^, i.e., V — 5L^, then 
the burning speed of the flame is deflned as = m/pfL^, 
where rh = e/q is the total fuel-consumption rate and pf 
is the fuel density. Then the condition e ^ e/ts can 
be rewritten as S" ^ CgE/qpf, where ts — S/cs, Cs is the 
sound speed, and E = s/V is the internal energy per unit 
volume. The flame here may be laminar, turbulent, or 
distributed, provided it has the required burning speed. 

In order to examine the physical meaning of this con- 
dition on S, assume an ideal gas equation of state, 
E — P/{-f — 1). At the start of the runaway, pressure 
is nearly constant across the flame. Then the prod- 
uct density is pp = PfTj/Tp = PfTf/{Tf + q/Cp) = 
P/{P/pf + q["f — l)/7), where Tp is the product tem- 
perature and Cp is the specific heat at constant pres- 
sure. For energetic reactive mixtures, the denominator 
P/ Pf + q{l— l)/7 can be approximated as q{"f — 1). Here 
q = 43.28i?To/M > Pq/po [i3| and at the onset of the 
runaway, P ~ 1.5Po and pf ^ po, giving the accuracy of 
this approximation »i6%. Thus, pp ~ P/qi"f — 1), and 
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where a = pf/pp is the fluid expansion factor. 

In the reference frame of a steady flame, PpUp = 
PfUf = PfS, where Uf and Up are the velocities of the 
fuel and product, respectively. Thus, eq. ([1} is equiv- 
alent to the statement that Up = Cs- If Cs is taken as 
the sound speed in the product, then the flame with the 
speed satisfying eq. (P) is a CJ deflagration [23 1. 

The speed of a CJ deflagration, Scj, is a theoretical 
maximum speed of a steady-state flame. The discussion 
above shows that such a flame generates enough energy 
on a sound-crossing time to raise its internal pressure 
and, thus, disrupt its steady-state structure. Real lami- 
nar flames, both chemical [23| and thermonuclear [l^, do 
not have burning speeds that approach Scj- Turbulent 
flames, however, can develop such high values of St- 

Unlike a laminar flame, in which the local sound speed 
increases smoothly from its value in the fuel, Csj, to that 
in the product, Cg^p, a turbulent flame effectively consists 
of two fluids with either Cgj or Cg^p- Figure O shows 
Scj based on both Csj and Cs,p. Dissipative heating of 
fuel by turbulence causes Csj and Cg^p to increase and 
a to decrease. Thus, the horizontal shaded gray areas 
show the range of values of Scj corresponding to fuel 
temperatures « 360 — 430 K. In particular, in Case 10, 
Tf Ri 360 K at 2Ted (lower bound of the shaded regions) 
and increases to «430 K by 6.5Ted (upper bound). 

Figure [2] shows that, upon first reaching the QSS, St 
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is close to, but still below, c^j/a, which prevents the 
onset of the runaway. During the time (2 — 6.5)red, 
turbulent heating of fuel increases Sl by a factor of 
~2, thus accelerating St above the critical value Cgj/a 
and allowing the runaway to begin. Figure [3Jc) shows 
that, at this point, the product velocity indeed becomes 
~ Csj- Furthermore, the growth rate of St increases 
significantly once St becomes >Cs_p/a, i.e., when Up be- 
comes supersonic relative to both sound speeds. Note 
also that the transition from a QSS to a detonation oc- 
curs on a sound-crossing time of the turbulent flame 
ts = 5t/cs.q ~ 27/iS ~ Ted, where (5t ~ 1 cm is the 
flame-brush width (Fig. [Sb) and Cs,o ~ 3.7 x 10'' cm/s. 

Figure [2] also shows St for turbulent H2-air flames for 
other values of Ui and I. In Cases 5-7, St remains well 
below Cs.f/a, and the flame evolves in the QSS, as de- 
scribed in [13, [13. This QSS was observed over signifi- 
cantly longer periods of time than shown in Fig. [2l e.g., 
WTed in Case 6. Cases 1-4 were similar and so are not 
shown. The runaway process was also observed in Cases 
8 and 9, in which, however, the flame accelerated quickly 
and left the domain before DDT could occur. Note that 
the overall growth rate of St in Cases 8 and 9 was lower 
than in Case 10 {red increases with decreasing Ui). 

To determine the dependence of the results on the reac- 
tion model, we carried out a similar simulation for a sto- 
ichiometric CH4-air mixture. In this case, 6l,o — 0.042 
cm is close to that in H2-air, but Sl,o = 38 cm/s is 
eight times lower llj. The CH4-air system also showed 
DDT, but at a higher turbulent intensity relative to Sl 
{Ui — 2.24 X 10"^ cm/s ~ 59Sl) and in a larger system 
{l = 0.31 cm, L = 1.328 cm) (Case 11, Figs. [T] and [2]). 
The overall evolution, however, was different from Case 
10. The time to DDT was « 2Ted, and the flame never 
developed a QSS. The flame accelerated significantly rel- 
ative to fuel, which required a longer domain to observe 
DDT, and, in contrast with Case 10, a strong well-defined 
global shock formed and ran ahead of the flame brush. 

The key aspect of the spontaneous DDT mechanism 
discussed here is that it does not place any specific con- 
straints on the equation of state, reaction model, or the 
flame properties. A decrease of fluid density with increas- 
ing temperature in an exothermic process means that, at 
a high but subsonic burning speed, the flow of products 
becomes supersonic relative to the flame, irrespective of 
how burning occurs. This ensures that the pressure wave 
remains coupled to the region in which the energy release 
occurs (note the location of peaks of P and E in Fig. [3)3 
and c). This is in contrast with the spontaneous reaction- 
wave model [l3| . which requires very specific hot-spot 
properties in order for the reaction wave and the pres- 
sure pulse it produces to remain properly coupled. 

Figure [1] suggests that there is both a minimal system 
size and a minimal relative turbulent intensity at which 
DDT is possible, and they appear to increase for reactive 
mixtures with slower laminar flames. Applying eq. ([1} to 



establish whether DDT can occur depends on our abil- 
ity to predict the turbulent flame speed for given Ui and 
I. This is particularly difficult in the high-speed regimes 
where spontaneous DDT is most likely to occur. Further 
studies using more detailed chemical kinetics models are 
required to establish the range of regimes in which DDT 
is to be expected for realistic reactive mixtures and to 
investigate the possibility of flame extinction in the pres- 
ence of high-intensity turbulence. 
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